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Abstract 

Homologous classes of Polycyclic Aromatic Hydrocarbons (PAHs) in their crys- 
talline state are among the most promising materials for organic opto-electronics. 
Following previous works on oligoacenes we present a systematic comparative study 
of the electronic, optical, and transport properties of oligoacenes, phenacenes, cir- 
cumacenes, and oligorylenes. Using density functional theory (DFT) and time- 
dependent DFT we computed: (i) electron affinities and first ionization energies; (ii) 
quasiparticle correction to the highest occupied molecular orbital (HOMO)-lowest 
unoccupied molecular orbital (LUMO) gap; (iii) molecular reorganization energies; 
(iv) electronic absorption spectra of neutral and ±1 charged systems. The excitonic 
effects are estimated by comparing the optical gap and the quasiparticle corrected 
HOMO-LUMO energy gap. For each molecular property computed, general trends 
as a function of molecular size and charge state are discussed. Overall, we find that 
circumacenes have the best transport properties, displaying a steeper decrease of 
the molecular reorganization energy at increasing sizes, while oligorylenes are much 
more efficient in absorbing low-energy photons in comparison to the other classes. 
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1 Introduction 



A large class of conjugated 7r-electron systems, polycyclic aromatic hydro- 
carbons (PAHs) are of fundamental importance in many research areas of 
chemistry as well as in materials science and astrophysics. PAHs are precur- 
sors of extented carbon networks |la|llbj and their carbon skeletons can be 
seen as small pieces of graphene |2af2b] . Some PAHs are of interest in en- 
vironmental chemistry due to their carcinogenicity and their ubiquity as air 
pollutants produced by the combustion of organic matter |3a|3b] . Neutral and 
charged PAHs up to large sizes are thought to be abundant in space and are 
thus seen as an intermediate stage between the gas and dust phases of the 
matter between stars |4af4bj . 

Homologous classes of PAHs in their crystalline state are among the most 
promising materials for organic electronics [5]. Oligoacenes and their deriva- 
tives, for example, are being increasingly used as active elements in a variety 
of opto-electronic devices such as organic thin-film field-effect transistors [6] , 
light-emitting diodes [7], photovoltaic cells [8J, and liquid crystals Organic 
electronics based on functionalized acenes and heteroacenes is presently a very 
active field of research jlOallObllOcllOdUDe] . 



The main features controlling organic semiconducting devices are exciton for- 
mation, migration, and dissociation, charge transport, charge collection at the 
electrodes, molecular packing in the bulk material, and absorption and emis- 
sion properties [14a|14bj . Since many of these properties pertain mostly to the 
molecular units, systematic quantum-chemical calculations on isolated PAHs 
can significantly contribute to our understanding of the electronic and optical 
properties of families with promising opto-electronic applications. 

As a part of a more extensive research on PAHs [15J, and following our 
previous work on oligoacenes |f 6a)16bj . we present in this paper a compre- 
hensive comparative theoretical study of four homologous classes of PAHs, 
namely oligoacenes, phenacenes, circumacenes, and oligorylenes, in their neu- 
tral, cationic, and anionic charge-states. This choice of families was moti- 
vated by the availability of reliable experimental data for the first members 
of each class. They include both catacondensed (oligoacenes and phenacenes) 
and pericondensed (circumacenes and oligorylenes) species, spanning a repre- 
sentative range of regular, relatively symmetric PAHs. Moreover circumacenes 
and oligorylenes converge, at the infinite limit, to a zigzag and arm-chair edged 
graphene nanoribbon, respectively, and as such are promising candidates for 
organic and molecular electronics [13] . We restricted ourselves to the first five 
members of each class. 

The geometries of the molecules considered are sketched in Fig. [TJ Oligoacenes 



2 



[5] and phenacenes [TT] consist of fused benzene rings joined in a linear and 
zig-zag arrangement, respectively. Oligorylenes [12] can be seen as a series of 
naphthalene molecules connected by two carbon-carbon bonds, and circum- 
acenes are obtained by addition of benzene units around the circumference 
of their oligoacenes counterparts. Note that, with the exception of circum- 
tetracene, circumpentacene, and hexarylene, all of the molecules considered 
have been synthesized. 

We used density functional theory (DFT) |17a|ll7b|ll7c|ll7d] to compute the 
electronic ground-state of the molecules considered, and time-dependent DFT 
(TD-DFT) |18a|ll8b|ll8c|ll8d] to obtain their electronic absorption spectra. To 
the best of our knowledge the complete photo-absorption spectra of higher 
neutral and charged phenacenes, circumacenes, and oligorylenes are reported 
here for the first time. The knowledge of this property, in particular, enabled a 
comparison of the absorption in the visible part of the spectrum, a quantity of 
interest for photovoltaic applications. Moreover, we used total-energy differ- 
ences between the self-consistent field calculations performed for the neutral 
and charged systems (±1) to evaluate: (i) the vertical and adiabatic elec- 
tron affinities and first ionization energies; (ii) the quasiparticle correction to 
the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular 
orbital (LUMO) gap; (iii) the molecular reorganization energies. The compari- 
son between the optical gap, i.e. the lowest singlet-singlet TD-DFT excitation 
energy, and the quasiparticle corrected HOMO-LUMO gap enabled a rough 
estimate of the excitonic effects occurring in these molecules. 

The paper is organized as follows. Section [2] reports some technical details 
of the calculations. The results we obtained are presented and discussed in 
Sect. |3j Our concluding remarks are reported in Sect. HI 




Fig. 1. Oriented geometries of the molecules considered. From top to bottom: 
oligoacenes (naphthalene, anthracene, tetracene, pentacene, hexacene), phenacenes 
(phenanthrene, chrysene, picene, fulminene, 7-phenacene) , circumacenes (coronene, 
ovalene, circumanthracene, circumtetracene, circumpentacene), and oligorylenes 
(perylene, terrylene, quaterrylene, pentarylene, hexarylene). 
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2 Computational details 

2.1 Geometry optimizations 

We used the Gaussian-based DFT module of the nwchem package [19] to 
compute the ground-state properties of each molecule in its —1, 0, +1 charge-state. 
Following our previous studies on linear acenes |16a|ll6b] , we used a basis-set 
with the smallest addition of diffuse functions, namely the 6-31+G* basis, a 
valence double zeta set augmented with d polarization functions and s and 
p diffuse functions for each carbon atom. Geometry optimizations were per- 
formed using tight convergence criteria, specified by maximum and root mean 
square gradient thresholds of 1.5-10 -5 and 1.0-10 -5 atomic units, respectively, 
and maximum and root mean square thresholds of the Cartesian step respec- 
tively of 6.0-10 -5 and 4.0-10" 5 atomic units. All neutral and singly-charged 
species were computed as singlet and doublet, respectively, which are known 
to be the ground states for this class of molecules [22b.23a,23cJ. 




Geometry optimization were performed using the hybrid B3LYP functional, 
a combination of the Becke's three parameter exchange functional |20j and 
the Lee- Yang-Parr gradient-corrected correlation functional [2T] . In compar- 
ison with other exchange-correlation functionals, the B3LYP has proven to 
give results markedly more accurate for the ground-state properties of PAHs 



a|22b] and has been used in a large number of studies |23a,23b 23c 23d 23o 231' 23g 231i 23i 23j ,23k 




Although the B3LYP functional has been developed for the electronic ground-state, 
it is also routinely employed to obtain electronic excited-state properties, and 
its application usually yields accurate results for low-lying valence-excited 
states of both closed-shell [24a|24bj and open-shell [25a }25b] species. It is 
known, however, that this functional show the wrong asymptotic behaviour, 
decaying faster than 1/R (i.e., exponentially) for large distances R from the 
nuclei. Among the well known and documented limitations of conventional 
functionals [18d] are: (i) the correct description of Rydberg |26af26bf26cj . 
doubly-excited |27af27bj . and charge-transfer [28a, 28b, 28c] excited-states, (ii) 
the failure for large, extended 7r-systems such as polyacetylene fragments and 
oligoporphyrins |29j, and (iii) the system-size-dependent errors found for the 
lowest short-polarized excited states of neutral oligoacenes [23f|. 

This latter error in particular, which is particularly relevant to the present 
study, was recently shown to be corrected by range-separated functionals 
|30j . While conventional hybrid functional such as B3LYP are defined with 
a costant fraction of Hartree-Fock exchange, the authors of Ref. [30] used 
range-separated functionals |31af31bf31cf31d] incorporating both a position-dependent 
fraction and an asymptotically correct contribution of Hartree-Fock exchange. 
The application of range-separated functionals to organic compounds show 
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that they provide, in general, a better accuracy than global hybrids [32a,32b,32c,32d ; 32e 



As to the molecular properties of interest in this work (electron affinities, ion- 
ization energies, quasiparticle corrected HOMO-LUMO gaps, exciton bind- 
ing energy, electronic excitations), we recently checked for oligoacenes both 
the sensitivity of our calculations to different exchange-correlation functionals 
|33a|3 3b] . and their reliability in comparison with the available experimental 
data. With the only exception of the ionization energy, which were found to be 
described better at the LDA level, we found that the hybrid B3LYP functional 
yields the overall best performance. To further asses this choice we performed 
some benchmark calculations using the range separated functional LC-BLYP 
|32e] , which was found to give the best results in describing the optoelectronic 
properties of oligoacenes [30]. Tabled] reports the comparison between B3LYP, 
LC-BLYP, and the available experimental data of oligoacenes. A similar com- 
parison for some selected properties for which we could find the corresponding 
experimental data is presented in Table [2] for the smallest species of the four 
PAH families considered, namely naphthalene, phenanthrene, coronene, and 
perylene. 

The results in Table [1] are consistent with those reported in Ref. [30J, the 
small differences being due to the larger basis set cc-PVTZ used in the latter 
work. As shown by the reported mean absolute errors (MAE) and root-mean- 
squared errors (RMS), the use of LC-BLYP is indeed (i) able to solve the 
system-size-dependent errors of B3LYP found for the lowest excited state 
corresponding to the HOMO-LUMO transition (p bands); (ii) gives ionization 
energies that are in much more agreement than B3LYP with experiments; and 
(iii) provides better quasiparticle corrected HOMO-LUMO gaps. At the same 
time, however, the B3LYP functional appears to outperform LC-BLYP with 
respect to the a and (3 bands, the electron affinity, and the exciton binding 
energy; the molecular reorganization energies are described with the same 
accuracy by the two methods. When considering different families as done in 
Table [2] the B3LYP results are always in better agreement with experiments, 
with the only exception of the ionization energy. Since definitive conclusions 
cannot be drawn from our benchmark calculations and a thorough comparison 
between the performances of global hybrids and range-separated functionals 
is beyond the scope of this work, we decided to extend our previous analyses 
for oligoacenes [16a||16bj . using the same level of theory B3LYP/6-31+G* for 
the other classes of PAHs considered in this work. 



2.2 Total-energy differences calculations 

From the structural relaxations performed for both neutral and charged sys- 
tems, we computed via total-energy differences the adiabatic electron affinities 
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Table 1 

Comparison between B3LYP, LC-BLYP, and the available experimental data of 
oligoacenes (number of benzene rings from 2 to 6). The three lowest electronic 
transitions (usually described by Clars's notation p, a, 0) are computed by means 
of linear response TD-DFT (see Sect. 12. 3p . Ionization energies, electron affinities 
ionization energy, quasi-particle corrected HOMO-LUMO gaps, exciton binding en- 
ergy, and molecular reorganization energy are computed via total energy differences 
(see Sect. 12. 2p . The experimental data are taken from the compilation in Ref. [16aJ 
with the exception of the molecular reorganization energies A taken from Ref. [14b . 
The last rows report the mean absolute errors (MAE) and root-mean-squared errors 
(RMS) for each quantity; all values are given in eV. 
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0.96 


5.85 
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0.133 
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and the adiabatic single ionization energies. At the optimized geometry of the 
neutral molecule we then evaluated the vertical electron affinity (EA V ) and 
the vertical first ionization energy (IE V ). This enabled the calculation of the 
quasiparticle-corrected HOMO-LUMO gap of the neutral systems; this quan- 
tity is usually referred to as the fundamental gap and is rigorously defined 
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Table 2 

Same as in Table Q] for the smallest species of each PAH family considered; the 
experimental data are taken from the references given in Tables OH 
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B3LYP 


4.36 
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8.27 
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3.63 




EXP 
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0.47 
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2.79 
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2.85 
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(0.13) 
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within the ASCF scheme |17bj as: 

QPgap — IE V — EA V = E N+1 + E N _! — 2E N , (1) 

E N being the total energy of the N-electron system. We used also the following 
approximate expression |17aj : 

QPgap = e N+l — e N' (2) 

where ej is the i th Kohn-Sham eigenvalue of the j-electron system. The results 
obtained using the above Eqs. (CQ) and ([2]) tend to coincide as the system gets 
larger and the orbitals more delocalized. 

Concerning transport properties, in organic molecular semiconductors con- 
ductivity is known to occur via a hopping mechanism in which charge carriers 
jump between adjacent molecules, usually under the effect of an external ap- 
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plied field. This process can be described by the semiclassical Marcus theory of 
electron transfer with overall rate constant proportional to exp [—\/(4KbT)] 
(K-q Boltzmann's constant, T temperature) |14a|ll4b| . The quantity A, the 
molecular reorganization energy or intramolecular coupling, affects critically 
the charge-transfer process: low values of the reorganization energy imply high 
transfer rates. A can be computed via the "four-point method" as the sum of 
two contributions [14a.l4bJ: 

X = \ 1 + \ 2 = (E1 -E+) + (E*-E). (3) 



In the equation above, E and E + are the total energies of the neutral and 
cation in their equilibrium structures, respectively, is the total energy of 
the cation in the neutral geometry, and E* is the total energy of the neutral in 
the cation geometry. While the first term Ai corresponds to the radical cation 
formation, the second one A2 accounts for the relaxation of the charged state. 



2.3 TD-DFT calculations 



To compute the excitation energies and electronic absorption spectra we used 
two different implementations of TD-DFT in the linear response regime, in 
conjunction with different representations of the wavef unctions: 

(1) the frequency-space implementation |18aj based on the linear combina- 
tion of localized orbitals, as given in the nwchem package [19] : 

(2) the real-time propagation scheme using a grid in real space |18bj . as 
implemented in the OCTOPUS computer program [35af35b| . 

We used the scheme ([T]) at the same level B3LYP/6-31+G* used to obtain 
the ground-state geometries and restricted ourselves to the first few singlet- 
singlet electronic transitions. In scheme ([2]) one obtains the whole absolute 
photo-absorption cross-section o~(E) from the dynamical polarizability a(E) 
through the equation: 

^j—QWi?)}, (4) 



where h is Planck's constant, is the imaginary part of the dy- 

namical polarizability, and c is the velocity of light in vacuum. The dipole 
strength-function S(E) is related to c(E) by the equation: 
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m e and e being respectively the mass and charge of the electron. S(E) has units 
of oscillator strength per unit energy and satisfies the Thomas-Reiche-Kuhn 

dipole sum-rule N e = J dE S(E), where N e is the total number of electrons. 

The real-time TD-DFT method in real space was proven to give good results 
for PAHs up to photon energies of about 30 eV [16a|ll6b| . For the calculations 
performed with OCTOPUS we used the same prescription as in Ref. [16aJ. 
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Fig. 2. Percentage variation of the rotational constants A (top panel) and B (bot- 
tom panel) relative to the neutral counterparts for anions (asterisks) and cations 
(crosses), as a function of molecular size for oligoacenes (black), phenacenes (green), 
circumacenes (blue) , and oligorylenes (red) . We omit the entries for naphthalene an- 
ion, which is known not to form a stable anion in the gas-phase. 
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3 Results and discussion 

3.1 Geometry optimization 

In our previous studies [ 16ajl6b] we found that the structural parameters ob- 
tained for ionized oligoacenes, ignoring any apparent symmetry and adopting 
the symmetry of the corresponding neutral molecule, are coincident within 
numerical errors. We therefore assumed in this work the same symmetry of 
the neutral molecule for each charged species considered, with the exception 
of coronene anion and cation, which are known to undergo Jahn-Teller dis- 
tortion from D 6 /j to D 2 h symmetry upon single ionization [31]. As previously 
done for oligoacenes [16aj . the structural variations occurring in the charged 
species compared to the corresponding neutral one are expressed in Fig. [2] in 
terms of the percentage variations of the rotational constants A and B. These 
latter quantities are proportional to the inverse of the principal momenta 
of inertia in the plane of the molecule. Since rotational constants are here 
used just as a compact way to represent structural changes between different 
charge states, they were all computed at the optimized geometries, disregard- 
ing vibration-rotation coupling effects. Evaluating the differences of rotational 
constants between the charged species and the corresponding neutral counter- 
parts, AA and AB, we thus computed the quantities A ro i = AA/A neutra i and 
B re i = AB/B neutral . Figure [2] shows that A rc i and B re i display well-defined 
trends function of molecular size. 

For oligoacenes and circumacenes the largest structural changes are observed 
for the anions, while the variations of cations are always small; in particular, 
the anions of both families appear to be primarily distorted along the long 
axis showing for all sizes a lengthening along it (B re i < 0) in comparison 
to the neutral molecules. The same holds true for both ±1 charge-states of 
phenacenes, the cations showing also a sensible shortening along the short-axis 
(A re i > 0). In the case of oligorylenes, instead, at increasing molecular size the 
largest structural changes are observed for the cations, which appear to be 
primarily distorted along the long axis showing a shortening along it (B rc i > 0). 



3.2 Ionization energies and electron affinities 

The adiabatic and vertical values of electron affinities and ionization energies 
as obtained via total energy differences at the B3LYP/6-31+G* level are given 
in Table [3J and displayed in Fig. [3J as a function of molecular size; for com- 
parison, we report also the available experimental data [36]. In particular, the 
differences between computed and measured first ionization energies, under- 
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Table 3 

Adiabatic and vertical values (within parentheses) of electron affinities and ioniza- 
tion energies of oligoacenes, phenacenes, circumacenes, and oligorylenes; all data, in 
eV, have been obtained via total energy differences at the B3LYP/6-31+G* level. 
The experimental data are taken from the NIST Chemistry WebBook |36j . 



Molecule 


Electron 


affinity 


Ionization energy 




Theory 


Experiment 


Theory 


Experiment 




Acenes (C 4n+ 2H 2ri+4 ) 






naphthalene 


-0.26(-0.38) 


-0.20±0.05 


7.80(7.89) 


8.144±0.001 


anthracene 


0.53(0.43) 


0.530±0.005 


7.02(7.09) 


7.439± 0.006 


tetracene 


1.08(1.00) 


1.067±0.043 


6.49(6.55) 


6.97±0.05 


pentacene 


1.48(1.41) 


1.392±0.043 


6.12(6.16) 


6.589±0.001 


hexacene 


1.78(1.72) 




5.83(5.87) 


6.36±0.02 




Phenacenes (C 4n+ 2H 2rH 


-4) 




phenanthrene 


-0.05(-0.21) 


0.01±0.04 a 


7.53(7.63) 


7.891±0.001 


chrysene 


0.29(0.19) 


0.397±0.008 


7.17(7.25) 


7.60±0.01 


picene 


0.40(0.10) 


0.542±0.008 


7.04(7.13) 


7.51±0.02 


fulminene 


0.57(0.48) 




6.88(6.95) 


7.36±0.02 


7-phenacene 


0.64(0.54) 




6.80(6.87) 


— 




Circumacenes (CsnH2 n +6) 




coronene 


0.47(0.38) 


0.470±0.090 


7.02(7.08) 


7.29±0.03 


ovalene 


1.17(1.11) 




6.36(6.41) 


6.71±0.01 


circumanthracene 


1.68(1.63) 




5.90(5.94) 




circumtetracene 


2.06(2.02) 




5.57(5.60) 




cir cump ent acene 


2.34(2.31) 




5.33(5.35) 






Rylenes (Ci nH 4n+4 ) 






perylene 


0.96(0.87) 


0.973±0.005 


6.57(6.64) 


6.960±0.001 


terrylene 


1.55(1.48) 




5.98(6.05) 


6.42±0.02 


quaterrylene 


1.91(1.85) 




5.62(5.68) 


6.11±0.02 


pentarylene 


2.17(2.11) 




5.37(5.43) 




hexarylene 


2.37(2.31) 




5.18(5.24) 





Taken from Ref. [37J. 
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estimated by theory in all cases, are of the same order of magnitude as in 
previous analyses [23a,23kJ: they increase systematically at increasing molec- 
ular size for each class (from ~4% for naphthalene to ~10% for hexacene, 
from ~5% for phenanthrene to ~7% for fulminene, from ~6% for perylene 
to ~8% for quaterrylene, and from ~4% for coronene to ~5% for ovalene). 
As already found in the literature [30J and confirmed in our benchmark cal- 
culations (see Sect. 12.11 Tables [Q12]), these errors in the ionization energies 
can be largely corrected by range-separated functionals. On the contrary, as 
compared to ionization energies, the computed electron affinities are found 
to be in better agreement with experiments in all cases, confirming the very 
good performances of B3LYP/6-31+G* to compute the electron affinities of 
PAHs [23b 23cf23ij . Oligoacenes display the largest variations for ionization 
energies and electron affinities as a function of molecular size. On the other 
hand, phenacenes show the smallest variations for the same properties. Cir- 
cumacenes and oligorylenes behave in the same way at increasing sizes, with 
ionization energies reaching the value of about 5.2 eV and electron affinities 



> 

o 



cn 
o 

LJ 



9 r 

8 r 
7 

6 \ 
5 L 



onization energies 



□ 

aO 4 



♦o 



♦ 

o 

o o 



o o 



A 

4 



Electron affinities 



O O 



O 

o 



r- 



10 20 30 40 50 60 70 

Number of carbon atoms 



Fig. 3. Computed adiabatic ionization energies and electron affinities as a function 
of size for oligoacenes (black triangles), phenacenes (green squares), circumacenes 
(blue circles), and oligorylenes (red diamonds). When available, the corresponding 
experimental values are represented by the filled symbols (cfr. Table [3]). 
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converging to about 2.4 eV. 



3.3 HOMO-LUMO gap and excitomc effects 

For each molecule considered, Table @] compares the HOMO-LUMO gap E™, 
obtained as difference of Kohn-Sham eigenvalues, the excitation energy of the 
HOMO-LUMO transition E™~ DFT as given by frequency-space TD-DFT, 
and the corresponding experimental value Eg*£. In Table |4] we compare also 
the quasiparticle-corrected HOMO-LUMO gap computed via Eqs. ([I]) and 
d2J) to the corresponding experimental value QP g ap = IE e xp-EA exp (cfr. Ta- 
ble E}. The theoretical exciton binding energy E bind is estimated through the 
difference QPg ap -E™ _DFT , and compared with its corresponding experimental 
value QP|ap - Egap- F° r eacn class considered all of these quantities are dis- 
played in Fig. H] and, as expected due to a reduction of quantum confinement 
effects, they all decrease as a function of molecular size. 

As shown in Table H] and Fig. H] in the case of oligoacenes the HOMO-LUMO 
gap obtained as difference between Kohn-Sham eigenvalues, with relative er- 
rors in the range 1-6%, gives a better description of the experimental optical 
gap as compared to TD-DFT. In the latter already mentioned be- 

fore, a system-size-dependent error is known to exist |23f] , with a relative 
error increasing from 2 to 20% when going from naphthalene to hexacene. For 
circumacenes and oligorylenes the two methods have approximately the same 
accuracy in predicting the optical gap, with relative errors in the range 1-8%. 
On the contrary, TD-DFT results for phenacenes agree better with experi- 
ment (relative error 1-4%) than the corresponding Kohn-Sham gap (relative 
error 8-10%). 

As previously discussed in Ref. [16a], the ASCF QP-corrected HOMO-LUMO 
gaps of neutral oligoacenes compare very well with the DFT-based tight-binding 
GW data but they are affected by similar errors as those found with 
TD-DFT. On the other hand, since these errors almost cancel each other in 
the evaluation of the exciton binding energy, we obtain an accurate estimate of 
Ebmd- To the best of our knowledge the QP-corrected HOMO-LUMO gap for 
phenacenes, circumacenes, and oligorylenes have never been reported before 
in the literature. By comparing these values with the corresponding optical 
gap predicted by TD-DFT we could estimate the excitonic effects occurring 
in these molecules. We found in all cases appreciable excitonic effects due to 
both quantum confinement and reduction of screening; Ebind decreases as the 
size of the molecule increases with similar slope for the different families and 
approaches the value of about 1.6-1.8 eV for the largest species considered. 
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Fig. 5. Molecular reorganization energies computed using Eq. [3] for oligoacenes 
(black triangles) , phenacenes (green squares) , circumacenes (blue circles) , and olig- 
orylenes (red diamonds), as a function of the inverse of the total number of carbon 
atoms Nc- 

3.4 Molecular reorganization energies 



The last row of Table H] reports for each molecule the molecular reorganiza- 
tion energies A obtained through Eq. [3j Our results for oligoacenes are in good 
agreement with previous DFT results |14a|ll4bj . Note that, while the perfor- 
mances of the B3LYP functional for the calculation of A in conjugated mate- 
rials can be improved by precisely calibrated DFT functionals [14bJ, its use 
in the present work is justified since it is expected to describe reasonably the 
differences between different families of PAHs. Due to strong electron-phonon 
coupling, the molecular reorganization energy A is known to be inversely pro- 
portional to the total number of carbon atoms Nc [40J. Figure confirms 
this general trend with the only exception of phenacenes, which display for A 
an oscillatory behaviour like for the other physical properties computed (see 
Figs.rjM]). For the different families we observe a markedly different behaviour: 
while the reorganization energy of oligorylenes decreases slowly as a function 
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of molecular size, circumacenes display a much more pronounced slope as com- 
pared to oligoacenes. These results are consistent with those previously found 
in Ref. [2b], and confirm that circumacenes are indeed good candidates for 
new emerging technologies. 
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Table 4 

Comparison between HOMO-LUMO gap, TD-DFT excitation energy of the HOMO-LUMO transition, quasi-particle corrected g; 
and molecular reorganization energy for the PAHs considered. All data have been obtained at the B3LYP/6-31+G* level of theory. 



lVTnlppiilp f formula ) 


E KS 

gap 


-pTD-DFT 
gap 


pexp 
J^gap 


OP 1 

^ gap 


OP 2 

Kl i ± gap 


np ex P 

^ r gap 


J- 1 bind 


-r^exp 
"^bind 




naphthalene (CioHg) 


4.74 


4.36 


4.45 a 


8.27 


8.12 


8.34 


3.91 


3.89 


0.178 


anthracene (C14H10) 


3.54 


3.21 


3.45 a 


6.66 


6.58 


6.91 


3.45 


3.46 


0.134 


tetracene (C18H12) 


2.74 


2.45 


2.72 a 


5.55 


5.50 


5.90 


3.10 


3.18 


0.109 


pentacene (C22H14) 


2.19 


1.91 


2.31 a 


4.75 


4.72 


5.20 


2.84 


2.89 


0.092 


hexacene (C26H16) 


1.78 


1.51 


1.90 a 


4.15 


4.13 




2.64 




0.079 


phenanthrene (C14H10) 


4.67 


4.19 


4.24 b 


7.84 


7.44 


7.88 


3.65 


3.64 


0.212 


chrysene (C18H12) 


4.21 


3.73 


3.89 b 


7.06 


7.00 


7.20 


3.33 


3.31 


0.161 


picene (C22H14) 


4.19 


3.70 


3.80 b 


7.03 


6.82 


6.97 


3.33 


3.17 


0.180 


fulminene (C26H16) 


4.00 


3.47 




6.47 


6.42 




3.00 




0.144 


7-phenacene (C30H18) 


4.00 


3.50 




6.33 


6.26 




2.83 




0.148 


coronene (C24H12) 


4.00 


3.42 


3.74 c 


6.70 


6.64 


6.82 


3.28 


3.08 


0.136 


ovalene (C32H14) 


2.91 


2.69 


2.88 c 


5.30 


5.28 




2.61 




0.095 


circumanthracene (C40H16) 


2.14 


2.03 




4.31 


4.30 




2.28 




0.073 


circumtetracene (C48H18) 


1.58 


1.54 




3.58 


3.58 




2.04 




0.056 


circumpentacene (C56H20) 


1.19 


1.18 




3.04 


3.04 




1.86 




0.042 


perylene (C20H12) 


2.98 


2.79 


2.82 d 


5.77 


5.71 


5.99 


2.98 


3.17 


0.142 


terrylene (C 30 Hi 6 ) 


2.22 


2.21 


2.22 d 


4.57 


4.53 




2.36 




0.128 


quaterrylene (C40H20) 


1.79 


1.83 


1.88 d 


3.83 


3.81 




2.00 




0.122 


pentarylene (C50H24) 


1.51 


1.54 


1.66 d 


3.32 


3.30 




1.78 




0.118 


hexarylene (C60H28) 


1.30 


1.36 




2.93 


2.92 




1.57 




0.117 



See compilation in Ref. [T6a]. 6 From Ref. [23e]. c From Ref. [38] ■ d From Ref. [12]. 



3.5 Absorption spectra 



Figures IMS display the absorption spectra of each of the molecule consid- 
ered, in its neutral and ±1 charge-state, as computed using the real-time 
real-space TD-DFT implementation [35a, 35b]. Note that the spectra of the 
higher neutral and singly-charged phenacenes, circumacenes, and oligorylenes, 
have never been reported before. This enabled a quantitative comparison of 
the absorption in the visible part of the spectrum for each class, which gives 
an estimate of its potential use as photo-active element in opto-electronic de- 
vices. To this end Fig. [9]shows the integrated values in the range 1.0-3.0 eV of 
the individual dipole strength-functions S(E) (see Eq. [5]) divided by the total 
number of carbon atoms in each molecule as a function of molecular size. 

As expected from previous analyses [16a|fl6 bJ. with the exception of the dif- 
ferent intensities which are related to the total number of electrons in the 
molecule, all of the spectra in Figs. |6]|8] display a similar broad excitation 
peaking at ~17-18 eV, which involves n — > a*, a — > n*, a — > a*, and Rydberg 
spectral transitions, and which is relatively insensitive to the charge-state of 
the molecule. Note that we already showed [16alll6b] our results to be in good 
agreement up to photon energies of about 30 eV with the experimental data 
obtained for a few neutral PAHs [38]. However, the use of a finite simulation 
box in our simulations does not give a satisfactory treatment of continuum 
effects in this spectral window and produces spurious structures as compared 
to experiments [16a|ll6bj . 

The main differences among families arise in the low-energy part of the spec- 
trum, from the near-IR to the near-UV up to about 8 eV, which involves 
electronic transitions of n — > 7r* character. This is also the spectral window 
which is more interesting from an applicative point of view. In particular, an 
inspection of Fig. [H] shows that each family displays a continuous redshift of 
the absorption at increasing sizes. This behaviour was expected as a conse- 
quence of the continuous decrease in bandgap already observed for each class 
(see Fig. H|). As to the charged species, as shown in Figs. [THE] we found that 
PAH radical cations and anions display intense optical transitions at lower 
energies than their parent molecule, and have very similar electronic spectra. 
These findings, already known for charged PAHs from extensive spectroscopic 
studies in frozen glassy organic solids [41a|l41bj . are in qualitative agreement 
with the particle-hole equivalence in the pairing theorem of Hiickel's theory. 

The behaviour described above at low energies translates for each molecule 
in each class into a systematic increase of the integrated strength-function in 
the range 1.0-3.0 eV when going from the neutral to the cation and anion 
(see Fig. [9]). This effect is more pronounced for oligoacenes and phenacenes 
and, as expected, the observed scatter between the different charge-states 
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Fig. 6. Electronic absorption spectra of the neutral molecules of the four PAH classes 
considered, as computed using the real-time real-space TD-DFT implementation 
of the OCTOPUS code. 



decreases with increasing molecular size. Interestingly, oligorylenes appear to 
be much more efficient in absorbing low-energy photons in comparison to 
the other classes; as an example, focusing on the neutral species with the 
same number of carbon atoms, the value of S(E)/Nq for quaterrylene exceeds 
by more than 80% the corresponding value for circumanthracene. Since red 
photons represent a significant portion of the solar energy, this result might 
have important technological implications. 
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Fig. 7. Same as Fig. [6] for the corresponding cations. 



4 Concluding remarks 



We presented a systematic theoretical study of the five smallest species of 
oligoacenes, phenacenes, circumacenes, and oligorylenes in their -1,0, and +1 
charge-states. From the ground-state structural relaxations performed at the 
B3LYP/6-31+G* level we computed the electron affinities, the first ionization 
energies, the molecular reorganization energies, and the quasiparticle correc- 
tion to the HOMO-LUMO energy gap. We found good agreement with the 
available experimental data as well as with previous theoretical results. Bench- 
mark calculations performed using the range-separated functional LC-BLYP 
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Fig. 8. Same as Fig. [6] for the corresponding anions. 

suggest that the agreement with the experiment can be further improved for 
some selected properties (e. g. the ionization energy); a systematic study will 
be the subject of future work. Concerning transport properties, circumacenes 
appear to have the most favorable properties for technological applications, 
displaying a steeper decrease of the molecular reorganization energy as the 
molecular size increases. We computed also the electronic absorption spectra 
using a compendium of the TD-DFT theoretical framework in both frequency 
space, to obtain the optical gap and thus estimate the exciton binding energy, 
and real-time real-space, to obtain in a single step the whole photo-absorption 
cross-section extending up to the far-UV. The results of this latter method, 
in particular, enabled a close comparison between the families considered with 
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Fig. 9. Comparison between the integrated values in the range 1.0-3.0 eV of the 
individual dipole strength-functions S(E) (see Eq.[5|) divided by the total number of 
carbon atoms Nc for the anions (asterisks), neutrals (triangles), cations (crosses) of 
oligoacenes (black), phenacenes (green), circumacenes (blue), and oligorylenes (red) 
as a function of molecular size. 

respect to their potential use as photoactive materials in opto-electronic de- 
vices. Overall, oligorylenes appear to be much more efficient in absorbing 
visible light and red photons in particular, which might have important tech- 
nological implications in the field of photovoltaics. 
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